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Support  Project  DMA  PE  6370 IB/3204/240. 
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INTRODUCTION 


Geoid  heights  derived  from  SEASAT  satellite  altimeter  observations  of  the 
ocean  surface  and  from  potential  coefficients  were  used  in  a  least-squares 
solution  for  determination  of  the  radius  and  origin  shifts  of  a  mean-earth 
ellipsoid.  This  report  gives  a  description  of  the  geoid  height  data  used  in 
the  solutions  for  ellipsoid  parameters,  results  of  the  solutions,  and  summary. 


SEASAT  ALTIMETER  OBSERVATIONS 


The  altimeter  observations  were  collected  by  the  SEASAT  satellite  that 
was  launched  26  June  1978  and  recorded  data  continuously  over  the  oceans.  The 
mission  ended  9  October  1978  after  approximately  1000  revs  of  altimeter  obser¬ 
vations  were  obtained.  These  observations  were  prefiltered  using  a  linear  fit 
and  a  fixed  sigma  as  criteria  for  discarding  possible  extraneous  points.  The 
prefiltered  observations  given  at  10  points  per  second  (pps)  were  averaged  to 
2  pps  and  reduced  to  best-estimated  along-track  altimeter  geoid  heights,  NA. 
Further  editing  included  examining  data  at  both  ends  of  time  continuous  seg¬ 
ments  and  eliminating  geoid  heights  with  corresponding  vertical  deflections 
outside  the  range  ±60  arc  sec.  The  reduction  of  the  altimeter  measurements  to 
geoid  heights  entailed  making  atmospheric  and  environmental  corrections,  com¬ 
bining  with  the  precise  orbital  ephemeris  (Colquitt  et  al. ,  1980)  and  applying 
a  Kalman  smoother  (West  et  al. ,  1977).  The  ephemeris  was  computed  using 
SEASAT  Doppler  observations.  Two  sets  of  altimeter  geoid  heights  (Table  1) 
were  selected  for  the  ellipsoid  parameter  solutions  from  the  available  revs  of 
data.  The  altimeter  geoid  heights  are  based  on  a  reference  ellipsoid  with  a 
radius  of  6378145  m  and  reciprocal  flattening  of  298.26.  A  plot  of  the 
subtracks  for  both  altimeter  data  sets  are  given  in  Figures  1  and  2. 

The  along-track  geoid  heights  were  averaged  to  give  a  mean  geoid  height 
NA  for  v  5°  x  5°  area.  The  earth's  surface  was  divided  into  2160  5°  x  5° 
areas  bounded  by  latitude  lines  ±75°  and  longitude  lines  0°  and  360°.  NA  is 
assumed  to  be  located  at  the  center  of  the  5°  x  5°  area. 


GRAVITY  DATA  FROM  SATELLITE  OBSERVATIONS 


Station  coordinates,  station  heights  above  mean  sea  level,  and  potential 
coefficients  derived  from  satellite  observations  were  used  in  the  determina¬ 
tion  of  the  separation  of  the  mean  sea  level  from  the  reference  ellipsoid,  the 
geoid  height  (Groeger,  1968).  This  separation  is  referred  to  here  as  the 
geoid  height  derived  from  potential  coefficients,  Ng.  a  subset  of  NWL9D 
station  coordinates  and  heights  (Table  2)  and  the  DoD  WGS72  (Seppelin,  1974) 
potential  coefficients  were  used  in  the  generation  of  Ns  at  1°  x  1°  grid 
points.  The  geoid  heights  are  based  on  a  reference  ellipsoid  with  a  radius 
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of  6378135  m  and  flattening  reciprocal  of  298.26.  Groeger  also  provides  an 
adapted  reference  ellipsoid  with  a  radius  of  6378138.58  a  and  flattening 
reciprocal  of  298.28  that  best  represents  the  equipotential  geodial  surface 


Table  1 . 

Revs  of 

SEASAT 

Altimeter 

Set  1 

Set  2 

Day 

Revs 

Day 

Rev  a 

207 

416-428 

211 

472-473 

475-477 

479-486 


208 

429-443 

209 

444-457 

212 

487-498 

500 

215 

530-541 

218 

573-576 

579-584 

216 

544-557 

217 

558-570 

572 

219 

587-597 

600 

220 

601-611 

613-614 

227 

228 

707-715 

716-727 

729 

221 

615-627 

222 

631-640 

643 

232 

773-784 

786 

224 

661-667 

669-670 

233 

787-801 

225 

673-686 

234 

802-812 

815 

235  816-819 

823-827 
829 


\ 


236 


830-832 


s»-  Si.- 


t  i  r  i  i  i  i  i  i  i - 1 - 1 - 1 

0  30  60  90  120  ISO  ISO  210  240  270  300  330  360 
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Figure  1.  142  Revs  of  SEASAT  Altimeter  Data  Selected  from  Days  207  to  225 


LONGITUDE 


Figure  2.  123  Revs  of  SEASAT  Altimeter  Data  Selected  from  Days  211  to  236 


Table  2.  Station  Positions  and  Heights  (NWL9D) 


Station 


Number 

Location 

Latitude1 

Longitude 

hmsz(m) 

he1 (m) 

8 

Brazil 

-23° 1 3 * 

314*  7* 

612.68 

605.27 

19 

McMurdo 

-77*50* 

166*40* 

38.2 

-17.84 

20 

Seychelles 

-4°40 ' 

55*28* 

592.5 

546.14 

21 

Uccles 

50047' 

4*21’ 

115.8 

148.  16 

22 

Philippines 

14°59' 

122®  4' 

11.5 

49.17 

23 

Guam 

13°26' 

144*38’ 

38.  18 

84.58 

24 

Samoa 

-14° 19 ' 

189*17* 

9.2 

39.07 

27 

Japan 

39°  8' 

141*  7* 

82.6 

114.12 

105 

Pretoria 

-25*56 ' 

28*20' 

1580.56 

1600.87 

107 

Virginia 

38*59* 

282*41* 

118.6 

76.86 

112 

Australia 

-34*40' 

138*39' 

34.44 

31.35 

113 

New  Mexico 

32*16* 

253*14* 

1205.84 

1172.78 

114 

Alaska 

61*17* 

210*10* 

67.6 

67.27 

116 

England 

51*11* 

358*37* 

78.93 

115.94 

118 

Greenland 

76*32* 

291*14' 

54.79 

59.14 

128 

Ottawa 

45*23* 

284®  4* 

86.08 

40.57 

192 

Texas 

30*23* 

262*16* 

178.5 

208.12 

311 

Maine 

44*24* 

291*59' 

25.88 

-11.68 

320 

Minnesota 

44*43* 

266*55* 

299.5 

259.75 

330 

California 

34*  6' 

240*56* 

461.56 

402.18 

340 

Hawaii 

21*31’ 

202®  0* 

401.19 

407.37 

641 

Florence 

43*48* 

11*13* 

100.3 

137.38 

9943 

Orroral 

-35*37* 

148*57' 

929.53 

942.56 

30078 

Massachusetts 

42*37' 

288*30' 

140.55 

104.15 

30120 

La  Paz 

-16*31’ 

291*49* 

4041.93 

4085.47 

30121 

Quito 

-  0*  5* 

281*34' 

2682.59 

2704.0 

30122 

Paraguay 

-25*18* 

302*23* 

177.16 

1893.86 

20284 

Sicily 

37*24* 

14*56’ 

20.83 

534.73 

30124 

Tehran 

35*44* 

51*23’ 

1421.22 

1415.76 

30126 

Zaire 

-4*22* 

15*15' 

453.88 

448.63 

30130 

Cyprus 

35*  0* 

33*43' 

101.9 

118.42 

30188 

Hawaii 

21*18* 

202*  0* 

4.01 

9.63 

30203 

Kenya 

-1*19* 

36*48’ 

1677.31 

1656.69 

30414 

Calgary 

50*52* 

245*42* 

1267.8 

1238.22 

30793 

Townsville 

-19*15* 

146*45’ 

6.41 

581.96 

30800 

Thailand 

13*47* 

100*35' 

15.34 

-22.32 

30967 

Bermuda 

32*19* 

295®  g> 

24.3 

-13.88 

1  Geodetic  coordinates  with  respect  to  an  ellipsoid  having  a 
radius  of  6378145  in  and  a  reciprocal  flattening  of  298.25 

2  Height  above  mean  sea  level 


Mean  geoid  heights,  Ns,  were  obtained  for  each  5*  x  5*  area  by  averaging 
the  Ns  values  at  the  1*  x  1*  grid  points  contained  within  the  area.  Ns  is 
referenced  to  location  X  at  the  center  of  each  5*  x  5*  area.  The  total 
surface  area  is  the  same  as  that  described  earlier  for  the  altimetric  geoid 
heights.  Geoid  heights  were  similarly  derived  from  the  GEM10B  (Lerch,  1978) 
potential  coefficients  using  the  same  station  coordinates  and  heights. 


SOLUTIONS  FOR  RADIUS  AND  ORIGIN  OF  A  MEAN-EARTH  ELLIPSOID 


Ellipsoid  parameters,  radius,  and  the  origin  of  a  mean  ellipsoid  best 
fitting  the  mean  geoid  height  differences  (AN)  obtained  from  altimetric  obser¬ 
vation  and  potential  coefficients  are  solved  for.  The  mean  geoid  heights  were 
derived  using  different  reference  ellipsoids.  If  it  is  assumed  that  the  ’tim- 
etric  geoid  from  the  ocean  areas  is  the  best  available,  then  in  order  tc  i.d 
the  equipotential  geoidal  surface  (as  defined  by  potential  coefficients)  at 
best  fits  the  difference  between  the  altimetric  and  potential  coeffieien 
geoids,  adjustments  are  made  to  the  geoid  derived  from  potential  coeffic  '.s 
such  that  the  mean  geoid  height  difference  is  zero.  Figure  3  shows  the 
sets  of  geoid  heights  and  the  associated  reference  ellipsoids.  NA  is  th 
mean  altimetric  geoid  height  located  at  the  center  of  the  ith  5*  x  5* 
square  and  based  on  reference  ellipsoid  EA,  Ng,  the  mean  geoid  height 
derived  from  potential  coefficients,  is  also  located  at  the  center  of  the 
ith  5*  x  5*  square  and  based  on  reference  ellipsoid  Es.  ga  and  Gs 
are  the  altimetric  and  potential  coefficients  geoids,  respectively. 


Figure  3.  Altimeter  and  Potential  Coefficients  Geoid  Heights 
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A  change  in  geoid  heights  from  different  reference  ellipsoids  and  a  shift  in 
the  ellipsoid  origin  may  be  expressed  as  given  in  Equation  1  (Heiskanen  and 
Moritz,  1967) 

„  „  8N  ,  3N  .  3N  ,  3N  . 

na  =  Ns  +  Ta-  Aa  +  tit  Ax  +  3y  Ay  +  rsr  Az  d  ) 

where  Aa,  Ax,  Ay,  and  Az  are  increments  in  the  ellipsoid  radius  and  origin 
solved  for  in  a  least-squares  fit  (see  Appendix  A)  to  the  geoid  height  differ¬ 
ences. 


The  adjustment  (Aa  +  (aA  -  as)]  derived  from  the  solution  is  applied 
to  the  radius,  "aj "  of  the  reference  ellipsoid  that  best  fits  the  geoid  from 
potential  coefficients.  Aa  is  the  discrepancy  between  the  altimetri-  geoid  and 
the  geoid  derived  from  potential  coefficients,  and  (aA  -  as)  is  the  differ¬ 
ence  between  the  radii  of  the  reference  ellipsoids  for  the  altimetric  and 
potential  coefficients  data,  respectively. 

The  mean  ellipsoid  with  radius  “a"  that  best  fits  the  equipotential 
geoidal  surface  derived  from  potential  coefficients  is  given  below  by  Equation 

2. 


a  =  a1  +  [Aa  +  (aR  —  ) ]  (2) 

The  results  of  the  solutions  are  shown  in  Table  3.  The  first  solution  was 
derived  from  geoid  heights  based  on  WGS72  potential  coefficients  and  altimetric 
data.  The  altimetric  geoid  heights  were  taken  from  142  revs  of  satellite  data 
selected  from  the  time  period  from  day  207  to  225.  Solution  1  is  weighted  by 
the  number  of  data  points  in  each  5°  x  5°  square.  Solution  2  is  a  repeat  of 
Solution  1  with  the  exception  that  geoid  heights  derived  from  WGS72  potential 
coefficients  were  replaced  by  those  obtained  from  GEM10B  potential  coef¬ 
ficients.  The  approximate  1  m  discrepancy  in  the  radius,  a,  of  the  mean-earth 
ellipsoid  as  shown  between  Solutions  1  and  2  is  due  to  the  difference  in 
potential  coefficients. 


Table  3.  Ellipsoid  Parameter  Solutions  in  Meters 


Solution 

Reference 

Geoid 

Number  of 
Points 

1 

WGS72 

992,602 

2 

GEM 1  OB 

992,602 

3 

WGS72 

884,453 

4 

GEM10B 

884,453 

Ax 

Ay 

Az 

a 

-0.09 

0.12 

-2.47 

6378134.30 

0.01 

0.19 

-2.55 

6378135.43 

-0.33 

-0.17 

-2.67 

6378134.17 

-0.22 

-0.07 

-2.65 

63781  35.29 
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Solutions  3  and  4  were  computed  from  the  second  set  of  altimetric  data 
(Table  1)  and  also  weighted  by  the  number  of  points  in  each  5°  x  5°  area. 


SUMMARY 


The  earth's  radii  derived  from  the  above  solutions  were  found  to  be  con¬ 
sistent  between  the  two  sets  of  altimetric  data  both  of  which  contained  mean 
geoid  heights  in  more  than  1750  5°  squares.  The  WGS72  and  GEM1 OB  potential 
coefficients  gave  a  difference  of  1.1  m  in  the  radii.  An  average  of  the 
results  from  the  two  sets  of  potential  coefficients  shows  the  reference 
ellipsoid  that  best  fits  the  equipotential  geoidal  surface  derived  from  poten¬ 
tial  coefficients  has  a  radius  of  6378134.9  m.  The  results  presented  in  this 
report  are  corrected  results  to  those  given  under  the  heading  "SEASAT-1  Alti¬ 
metric  Value  for  Earth's  Semimajor  Axis  (West)"  by  Anderle  (1980). 

The  earth's  radius  has  been  derived  by  others  (Anderle,  1980)  from  sev¬ 
eral  data  sets  with  results  that  range  from  6378139  ±1m  derived  from  laser 
station  positions  to  6378134.5  m  derived  from  GEOS-3  satellite  altimetry  data. 
The  radius  of  the  mean  earth  reference  ellipsoid  obtained  here  falls  within 
the  above  range.  The  results  of  this  study  contributed  to  the  value  of 
6378136  m  recommended  as  a  standard  value  by  Anderle  (1980)  to  the  Interna¬ 
tional  Union  of  Geodesy  and  Geophysics  in  December  1979.  At  this  meeting  a 
new  Geodetic  Reference  System  1980  with  the  radius  of  the  earth,  a  =  6378137  m, 
was  adopted. 

The  origin  shifts  shown  in  Table  3  are  less  than  1  m  for  the  x  and  y  com¬ 
ponents.  An  average  of  the  shifts  in  the  z  components  from  both  WGS72  and 
GEM10B  potential  coefficients  is  -2.5  m.  If  the  shift  is  directly  due  to  an 
error  in  the  origin  of  the  Doppler  Tracking  Station  Coordinate  system,  this 
shift  means  that  a  value  of  2.5  m  should  be  added  to  the  z  component  of  the 
station  coordinates.  In  a  report  by  Marsh,  (1980)  the  existence  of  a  4  m 
difference  between  Goddard  Space  Flight  Center  and  the  Naval  Surface  Weapons 
Center,  z  components  of  tracking  station  coordinates  is  presented.  The  -2.5  m 
shift  determined  in  this  paper  is  a  possible  explanation  for  part  of  that 
difference. 
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APPENDIX  A 

LEAST-SQUARES  FITTING  PROCEDURES 


A- 1 


LEAST  SQUARE  FITTING  PROCEDURE 


V 


The  given  data  points,  ANif  are  defined  as 
AN^Xj.)  »  NAi(Xi)  -  Nsi(Xi) 

the  difference  between  the  altimeter  geoid  heights,  and  geoid  heights  derived 
from  potential  coefficients  at  location  Xi  ■  Xi  ($ir  Xj)  for  i  -  1,2,...n 
where  n  is  the  total  number  of  observations.  These  points,  ANif  are  fitted 
to  the  geoid  height  difference  function,  AN*,  given  below. 

AN*  (Pj  ,x^)  ■  Aa  +  cosi^  cosXjAx  +  cos^sinX^y  +  sin^Az  (A-1) 

where  Pj (j-1 ,2,3,4)  are  the  ellipsoidal  parameters;  Aa,  Ax,  Ay,  and  Az, 
respectively.  These  are  increments  to  the  radius  and  origin  components  of  the 
reference  ellipsoid  for  the  geoid  derived  from  potential  coefficients. 

The  purpose  here  is  to  obtain  the  best  estimates  for  the  unknown  ellipsoid 
parameters  stated  above. 

ANa  «  AN* (Pj ,Xi)  +  eA  (A-2) 

The  "Principle  of  Least  Squares"  criterion  for  selecting  a  particular  solution 
to  Equation  (A-2)  is  to  determine  the  best  (Pj),  which  minimizes  the  sum  of 
errors  squared  after  the  solution 


. — 

FREC&UUQ  UtfNMflf 

A-3  L 


S  =  Z  VLel  =  ANt  -  ANMPj^)]  2 

wi  is  the  weight  matrix  as  defined  by  the  number  of  geoid  height  values  used 
in  computing  the  mean  geoid  height  for  the  ith  5°  x  5°  area. 

At  this  point  it  is  convenient  to  restate  the  problem  in  matrix  notation. 
Given  the  equation 

AN  =  AAP  +  E 

with  weight  matrix  W  and  A  =  .  AN  is  the  observed  data,  AN*  is  the 

computed  data,  and  AP  represents  the  ellipsoid  parameters,  Pj.  Find  the 
particular  solution  that  minimizes  the  sum  of  errors  squared  after  the 
solution, 

S  =  e’Vie  =  (AN  -  AAP)TW(AN  -  AAP) 
the  solution1  is  given  by 

AP  =  (a1Wa)“1a'iwan 

where 

T 

B  =  A  WA 
C  =  A^AN 
AP  *  B_1C 


1  John  L,  Junkins,  "An  Introduction  to  Optimal  Estimation  of  Dynamical  Systems" 
(Netherlands,  Sijthoff  and  Noordhoff  International  Publishers,  B.  V.  Alphen 
aan  den  Rijn,  1978). 


In  expanded  form,  the  matrices  are  given  as  follows: 


i 


*13  “  i£1co8^i3i„xiWi  b24  -  j^8 in^cos^cosJ^ 


S14  =  b33  =  jj  co^^sin**^ 

n  is  the  number  of  S®  v  «;•  ,  . 

areas  and  j  is  the  number  of  parameters. 
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